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ABSTRACT 

A STUDY OF BOILING WATER FLOW REGIMES AT LOW PRESSURE 
Mario P. Fiori 



Submitted to the Department of Nuclear Engineering in partial 
fulfillment of the requirement for the joint degrees of 
Nuclear Engineer and Master of Science in Mechanical Engineer- 
ing . 



A comprehensive experimental program to examine flow 
regimes at pressures below 100 psia for boiling of water in 
tubes was carried out . 

An electrical probe, which measures the resistance of 
the fluid between the centerline of the flow and the tube 
wall, was used to identify the various flow regimes o This 
probe proved to be an ideal detection device, because of its 
simplicity, reproducibility, and accurate representation of 
the flow pattern within the heated test section. 

The major flow regimes observed were bubbly, slug and 
annular flow. Under certain conditions at high flow rates, 
a wispy-annular flow pattern was observed. The effects of 
mass velocity (0.2 x 10° - 2.4 x 10^ lbm/hr-ft2), inlet temp- 
erature (100, 150 , 200^F), exit pressure (30, 100 psia), 
quality (x = -10 - +7 percent)., purity (9, 40 PPM NaCl; 1-3 
megohm-cm), length (L/D=<30, 60, 90)^ diameter (0,094, 0.242 
in.), and orientation (vertical and horizontal) on the flow 
regimes were studied. Flow regime maps on coordinates of 
mass velocity and quality are presented for these conditions. 

Bubbly and slug flow occurred primarily in the sub° 
cooled region, while fully developed annular flow was reached 
at equilibrium qualities between 2 and 4 percents The tran- 
sitions between the different flows were shifted to regions 
of Increased subcooling when velocity, pressure, and heat 
flux increased, and when inlet temperature decreased. Purity 
and geometry had little affect on the flow regime boundaries. 

The shifting of the transitions is related to the agglom- 
eration point, which is that point at which the bubbles so 
coalesce that slug flow is first observed. The agglomeration 
point depends on the point of Incipient boiling, the number 
of bubbles in the flow, and the number of collisions per 
bubble. These latter quantities in turn depend on velocity, 
temperature, pressure, and heat flux. 

The flow regime information obtained in this study should 
be of value in correlating and interpreting low pressure heat- 
transfer data. The flow regime data were found to be useful 
in explaining the effect of inlet temperature on burnout heat 
flux . 

Thesis Supervisor: Arthur E. Bergles 

Title: Assistant Professor of Mechanical Engineering 
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NOMENCLATURE 

A = heat transfer area (ft^) 

Cp = specific heat (Btu/lbtn^P) 

D = diameter (in.) 

p 

G = mass velocity (ibra/hr-ft ) 

h = heat transfer coefficient (Btu/hr-ft^°F) 

hg = enthalpy at exit (Btu/lbm) 

h^ = enthalpy of saturated liquid (Btu/lbm) 

hj^g = latent heat (Btu/lbm) 

h^^ = enthalpy at entrance to test section (Btu/lbm) 

I = shunt voltage (mv) 

k = thermal conductivity (Btu/hr-ft^F) 

L = length (in., ft) 

F = pressure (psia) 

= pressure at exit of test section (psia) 
q = rate of heat flow (Btu/hr) 

q/A = heat flux (Btu/hr -ft^) 

(q/A)^ = heat flux at incipient boiling (Btu/hr-ft^) 
(q/A)cr= critical heat flux (Btu/hr-ft^) 

= oscilloscope internal Impedance (ohms) 

R^ = resistance between probe and tube wall (ohms) 

R^ = variable resistance (ohms) 

g = water temperature at inlet or exit of tube (‘^F) 
= water bulk temperature (^F) 

= wall temperature (°F) 

= saturation temperature (°F) 

= voltage across test section (volts) 
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Vg == battery voltage (volts) 

Vj = voltage across variable resistor (volts) 
w = mass flow rate (ibra/hr) 

X = quality 

|i = viscosity (ibm/hr-ft) 

Plovj Regime Designations 
B = # = bubbly 

BTS = ® = bubbly- to- slug transition 

S r= o = slug 

STA = a = first indication of slug-to-annular transition 
STA = A slug-to-annular transition 

A - ■ = annular 

WA = jjJ) = wispy annular 



Chapter I 
INTRODUCTION 

1 .1 Background of Problem 

In recent years the interest in two-phase, gas-liquid 
flow has dramatically expanded as is indicated by the tre- 
mendous increase in the number of experiments and published 
papers in this field. If one examines the literature, one 
will find that a relatively large number of models have 
been proposed for analyzing and correlating experimental 
pressure-drop, heat-transfer, critical heat flux, and void- 
fraction data. These models frequently depend on the flow 
regimes (or flow patterns) which are geometric descriptions 

of the distribution of the liquid and gas phases within a 
/ 1 \ ^ 

channel. Gouse' ' considers knowledge of the flow regime 

to be as important as knowing whether the flow is laminar 

or turbulent in single-phase flow. 

/ 2 ^ 

Baker' ' showed that it was not practical to disregard 

the flow regimes in devising a two-phase pressure-drop 

( 3 ) 

correlation. Martinelll ' s' ' assumptions in his pressure- 
drop correlations tend to limit them to annular flow, 

Lavin and Young^^^ propose boiling heat-transfer correla- 
tions for refrigerants, which are dependent on the flow 
regimes. Numerous models for the critical heat flux condi- 
tion are based on the postulate of spray annular flow; 
however, there are some analyses which assume that fog flow 
is present. 



*Numbers in parenthesis refer to References on p. 80. 
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It is apparent, then, that it is important to know 
the flow regime if one is to utilize an appropriate model; 
however, direct relationships between flow regimes and 
observables, i.e., pressure drop, heat-transfer coeffi- 
cient, critical heat flux, and void fraction have not been 
clearly established Itistill appears, though, that 
flow regime observations will be helpful in reducing the 
complications of correlating the quantities of basic 
Interest . 

A survey of the literature Indicates that -flow regime 
studies have almost always been conducted under adiabatic 
conditions. However, exploratory studies have shown that 
flow patterns are considerably altered when there is heat 
addition, particularly in low pressure systems. 
Comprehensive flow regime maps have been prepared at 
Dynatech^'^^ and Harwell^^^ for high pressure boiling 
water systems. No similar study has been conducted at low 
pressures (below 100 psla). 

Low pressure systems are receiving increased atten- 
tion, particularly in the nuclear reactor and desalinlza- 
tion programs. In a recent study of critical heat fluxes 
at low pressure, Lopina^^^^ observed a pronounced effect 
of inlet temperature. For constant geometry, flow rate, 
and pressure level, the critical heat flux Increased as 
the inlet temperature increased. This effect, which is 
contrary to what is observed at high pressures, appears to 
be related to the flow regimes. The phenomenon of 
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upstream boiling burnout in uniformly heated channels also 
appears to be related to flow regimes, in particular, the 
slug flow regime, 

1 .2 Flow Regimes 

There have been various descriptions of flow regimes 
by different Investigators, These classifications have 
ofter been quite arbitrary, and are also complicated by the 

fact that vertical and horizontal flows cause different 
( 2 ) 

regimes. Baker^ ^ described no less than seven different 
regimes in adiabatic two-phase, gas-liquid flow in horizon- 
tal pipes. Vohr^^^^ and Kepple and Tung^^^^ have presented 
complete bibliographies on studies of adiabatic flow 
patterns , 

The flow regimes which may be encountered in flows in 
tubes with vertical upflow can be described in the order of 
Increasing gas or vapor content. If very little gas or 
vapor is Involved, the flow will be characterized by isola- 
ted vapor or gas bubbles dispersed rather uniformly in the 
continuous liquid phase. This is kn,own as ’'bubbly flow”. 
More bubbles appear with an increase of vapor or gas until 
the bubbles begin to touch one another and agglomerate. 

As more bubbles agglomerate, the gas begins to flow in 
long cylindrical, bullet-shaped bubbles, Between the 
bubbles are slugs of liquid and hence "slug flow". 

Further Increase of the gas content destroys the stability 
of the bubbles, producing "churn flow" or "semi-annular 
flow" in which slugs of liquid separate adjacent gas 



bubbles only temporarily. As the gas content continues to 
increase, the liquid flows on the wall as an annular film 
and the gas or vapor flows in the center of the tube. This 
geometry is called "annular flow". "Fog flow" or "dispersed 
flow" then follows where the gas or vapor phase occupies 
the entire tube but small drops of liquid are dispersed 
within it. This flow is also called "mist flow". 

Other flow patterns which have been encountered are 
"froth flow", "spray annular flow" and "wispy annular flow". 
Froth flow is quite similar to bubbly flow except the 
bubbles are so numerous that the flow has a frothy appear- 
ance . Spray-annular flow is a combination of mist and annu- 
lar flow. Some of the liquid flows in an annular film and 
some flows as small droplets in the vapor core. Wispy annu- 
lar is similar to spray annular except the flow in the core 
consists of agglomerated liquid instead of small droplets. 

There are very similar flow descriptions for horizontal 
flow. However, stratification is sometimes encountered in 
horizontal flow and several new regimes are introduced. 

"Plug flow" is similar to slug flow except that the bubbles 
move along the upper part of the channel or pipe . In 
"stratified flow" the liquid flows at the bottom, and gas or 
vapor at the top of the pipe. 

These regimes appear to be quite well defined. However, 
for practical reasons it is desirable to minimize the 
number of regimes. Gouse^^^ suggests that there are only 
four distinct flow patterns for vertical upflow: bubbly. 
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sliig^ annular, and mist; -while all the others which have 
been observed are merely transitions from one to another. 

The above flows were described for adiabatic condi- 
tions. In heated systems such clean geometric descriptions 
of the flow are not expected because of several character- 
istics not encountered under adiabatic conditions. Some of 
these characteristics are: 

1) Vapor is continually formed along the tube. This 
creates a volume flow rate which continuously increases 
along the tube; thus the flow regimes would be expected to 
be different in different sections of thetube . 

2) Bubbles may be generated and disturb the flow 
pattern. 

In any case, the adiabatic flow regimes are a basis for 
Interpreting and classifying the flow regimes with heat 
addition. The classification and both qualitative and semi- 
quantitative description of the flow patterns observed in 
this study are discussed in Section 2.3.3. 

1 .3 Flow Regime Detection 

Numerous methods of studying flow regimes in heated 
systems have been developed. High-speed still and movie 
photography has been most popular. More exotic methods 
which have been used Include isokinetic flow sampling, 
gamma or beta ray attenuation, and hot-wire anemometry. 

Each of these methods has drawbacks. Either there is too 
much data to analyze, as in the case of movies, or the 
methods are very complicated and expensive . The electri- 
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cal resistance probe which was 'developed by Solomon^ 
specifically for flow regime detection, is an inexpensive, 
simple, and reliable device. 

This probe measures the resistance between the center- 
line of the fluid flow and the tube wall. A bridge of 
liquid across the tube will cause a closed circuit, and 
gas or steam across the tube will result in an open circuit. 
Nassos^^^^ and Haberstroh^^^^ have successfully used this 
type of probe for adiabatic flow. Griffith^ adapted 
this probe for a high pressure boiling water system and 
delineated the slug-to -annular transition. Dynatech^"^^ 
has recently been using this technique to map boiling water 
flow regimes at high pressures (500-1000 psla)„ 

1 .4 Scope of Research 

A comprehensive study to examine boiling water flow 
regimes in tubes of small diameter at pressures below 100 
psia was undertaken. A resistance probe described in 
detail in Section 2.2 was chosen as the primary means of 
flow regime detection; however, a visual section was also 
employed. The experimental program consisted of an inves- 
tigation of the effects of velocity, inlet temperature, 
exit pressure, quality, (heat flux), purity, length, diam- 
eter, and orientation on the flow regimes. The first 
phase of the study was devoted to probe development and 
interpretation. 
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Chapter II 

EXPERIMENTAL PROGRAM 
2 .1 Description of Apparatus 
2,1.1 Hydraulic System 

The experimental facility used was the low-pressure 
test loop located In the MIT Heat Transfer Laboratory. The 
basic apparatus was designed and constructed In 1961.^^^^ 

A schematic of the loop is presented In Fig, 1, The 
pipings and fittings are made of brass and stainless steel 
for corrosion resistance. Rayon reinforced rubber hose is 
used where flexible connections are required. Distilled 
water is circulated by a bronze, two-stage, regenerative 
pump providing a discharge pressure of 260 psig at 3.6 gpm. 
The pump is driven through a flexible coupling by a 3-hp 
Allls-Chalmers induction motor. A Pulflo filter is in- 
stalled at the pump inlet. Pressure fluctuations at the 
outlet of the pump are damped out by means of a 2.5-gal 
Greer accumulator charged with nitrogen to an Initial 
pressure of 40 psig. This accumulator contains a flexible 
bladder-type separator which prevents the nitrogen from 
being absorbed by the system water. After the accumulator, 
the flow splits into the by-pass line and the test-section 
line . 

In the test-section line, fluid flows through a 
Flscher-Porter flowrator followed by a preheater, thence 
through a Hoke metering valve and the test section, after 
which it merges with fluid from the by-pass line. The 
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flovi then goes through the heat exchanger and returns to 
the pump. The preheater consists of four Chromalox heaters 
of approximately 6 kw each. Three of these are controlled 
simply with "off-on" switches while the fourth can provide 
a continuous range from 0 to 6 kw by means of a bank of 
two variacs mounted on the test bench. Quick-action James- 
bury ball valves are installed before the inlet to the 
flowrator and after the exit from the test section. This 
permits quick isolation of the test section to minimize 
fluid loss when conducting burnout tests. The exit valve 
is also used to adjust the test-section pressure. 

Flow through the by-pass line is controlled by a ball 
valve on each side of which there is a 300-pslg pressure 
gage. Pump operating pressure, and hence the pressure up- 
stream of the test section, is controlled by this valve . 

The heat exchanger is a counterflow type with system 
water flowing in the inner tube and city water in the 
outer annulus. Except at very high power levels, the heat 
exchanger maintained a constant pump inlet temperature. A 
Fulflo filter is installed on the city water line to 
reduce scale formation in the exchanger. 

The distilled water was de-ionized continuously, 
except during contamination studies, by passing a portion 
of the flow through four mixed-bed resin demlnerallzer 
cartridges Installed in a Barnstead- "Bantam" demineralizer 
unit. A 4.7-gal degassing tank is provided with five elec- 
trical heaters (3-200 vac and 2-110 vac). This tank also 
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serves as a surge tank. A 15~gal stainless-steel storage 
tank for filling the system is mounted directly above the 
degassing tank and can be filled with distilled water from 
standard 5-gal bottles with a small Hypro pump. 

2.1.2 Power Supply 

Test-section electrical power is provided by two 
36-k\j dc generators connected in series. Each generator is 
rated at 12 volts and 3000 amperes. The power control con- 
sole permits coarse or fine control from 0 to 24 volts, A 
water-cooled shunt Installed in parallel with the test sec- 
tion protects the generators against the shock of the sudden 
open circuit which occurs at burnout , Power is transmitted 
from the main bus to the test section by water-cooled power 
leads. In the present tests, the downstream or exit power 
lead was at ground potential. Rubber hose connected both 
inlet and exit chamber plenums to the main loop in order 
to electrically Isolate the test section. The exit of the 
test section was then separately grounded (see schematic 
of test section. Fig. 2). 

2.1.3 Ins t rume nt a 1 1 on 

All temperatures were measured by copper-constantan 
thermocouples made from 30-gauge Leeds and Northrup duplex 
wire. The test-section inlet temperature was measured by 
a thermocouple directly in the fluid stream, upstream of 
the Hoke metering valve. The thermocouple was introduced 
at this location through a Conax fitting equipped with a 



lava sealant . 
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All pressures were read on Bourdon-type gages located 
as shown in Fig. 1. The test-section inlet and exit pressures 
were measured with Helicoid 8-1/2 in. gages of 200 psig and 
100 psig, respectively. Both are specified to an accuracy 
of ^0.25^ of full scale. 

A variety of metering tubes and floats, which could be 
Installed interchangeably in the basic Fischer-Porter flow- 
meter housing, provided measurement of the test-section flow 
from 1.5 to 4000 Ibm/hr . 

The voltage drop across the test section was read 
directly on a Weston multiple -range dc voltmeter with a spec- 
ified accuracy of + 1/2 percent* The current flow was 
determined by using the Minneapolls-Honeywell, Brown recor- 
der to measure the voltage drop across a calibrated shunt 
(60.17 amp/mv) in series with the test section. 

The details of the resistance probe are discussed in 
Section 2.2. Polaroid cameras were used for the photo- 
graphic study. The camera settings and procedures are dis- 
cussed in Appendix B. 

2.1.4 Operating Procedures 

After the test section was installed, the loop and de- 
gassing tank were filled with distilled water from the 
supply tank. The water in the degassing tank was then 
brought to a boil while the loop water was circulated with 
the heat-exchanger coolant off. Degassing was accomplished 
by by-passing a portion of the loop water into the top of 
the vigorously boiling degassing tank. This was continued 
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until the temperature of the loop rose to approximately 
180^P. A standard Winkler analysis described in Ref. 17 
Indicated that this method of degassing reduced the air 
content to less than 0.1 cc air/liter. After completing 
degassing, the loop water was adjusted to the desired in- 
let temperature and the desired pressures were set. A 
pressure of from 210 psig to 260 psig was maintained at all 
times on the upstream side of the Hoke metering valve to 
prevent system-induced instabilities. The generators were 
then started and allowed to warm up. 

Nine burnout tests were made in order to have reliable 
critical heat flux data for conditions similar to the planned 
flow regime studies. This knowledge aided in avoiding burn- 
out of the probe test sections. Power was applied to the test 
section in small increments while maintaining a constant flow 
rate and inlet temperature. The test-section flow rate, 
shunt voltage (test-section current), inlet temperature, and 
exit pressure were recorded for each voltage setting. The 
critical or burnout heat flux was defined to be that heat flux 
which caused physical failure of the test section. After 
commencing flow regime studies, it was seen that the area of 
interest was usually far removed from (q/A)^^. Consequently 
no further burnout studies were made. Appendix A contains a 
tabulation of the burnout points. 

The same data were recorded for the flow regime runs. 

A run consisted of taking data at a particular mass flow 
rate for a given set of conditions, i.e. D, L/D, T^, . 
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As the power was increased^ various responses on the oscil- 
loscope indicated the type of flow (see Probe Interpretation, 
Section 2,3). Water purity was checked by means of a 
Barnstead Purity Meter (0.1 - 18 megohm-cm) or an Industrial 
Instruments Solu-Bridge Meter (0 - 30 PPM NaCl). 

2.2 Resistance Probe and Test Section 

2.2.1 Principle of Operation 

The resistance probe was used to determine the flow 
regime at the exit of the test section. The probe is a wire 
placed in the center of the tube as shown in Fig. 2. This 
wire is covered by teflon sleeving except for the tip of the 
wire. The circuit diagram for the probe is shown in Fig. 3a 
When there is a bridge of water across the tube at the test- 
section exit, there is essentially a closed circuit between 
the probe and the tube wall. This results in a maximum volt 
age across resistor which is measured by an oscillo- 
scope , On the other hand, steam passing between the probe 
and the wall resembles an open circuit giving a voltage 
reading across R^ near zero, 

2.2.2 Probe Integration with Test Section 

Two types of test sections were used for this study. 

The first type was used for preliminary burnout studies. 

The second type incorporated the resistance probe and a 
sight section. The sight section was not always used 
because the rubber gasket would not hold for exit pressures 
above 40 psig. The heated sections were made of standard 
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304 stainless tubes. The burnout test sections and the first 
probe test section had an 0.242-ln. l.d. and 0.3125-ln. o.d. 
while the second probe test section had an l.d. of 0.094 in. 
and o.d. of 0.120 In. The details of the two test sections 
are shown In Fig. 2. 

In the burnout test section, 3/4-ln. brass bushings 
were silver-soldered to the stainless steel tube to form 
the power connections upstream and downstream. A pressure 
tap was installed in each exit bushing with a No. 80, 
0.0135-ln. diameter hole drilled through the tube wall. 

The pressure tap was placed here rather than in the exit 
plenum chamber because of the choking phenomena whereby the 
plenum pressure is lower than the exit pressure . 

The probe test section had pressure taps located at 
both the inlet and exit. The exit section was modified to 
include the probe and pyrex sight section. The inside 
diameter of the glass tubing was equal to the outside diam- 
eter of the test section. This caused a slight expansion 
at the entrance to the visual section. 

The brass supporting structure for the probe, pressure 
tap, and sight section also served as the exit power connec- 
tion. The pressure tap is similar to the previous one 
described . 

The probe was fabricated from 0.035~ln. diameter temp- 
ered stainless-steel spring wire for the 0.242-ln. diameter 
tube and from 0.015-ln. wire for the 0.094-ln. diameter 
tube . Teflon spaghetti was slipped over this wire and 
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stretched until the teflon set snugly on the wire . The 
final probe diameters were O.O 6 O in. and 0.031 in. Even the 
small diameter probe had adequate stiffness to prevent vi- 
bration which could create spurious signals. 

The straight probe was pushed into the test section 
from the packing gland side until it touched the opposite 
wall. The tip was then carefully bent and the teflon was 
forced around the wire until only a 1/4 in. of probe tip was 
bare . 

The teflon acts as an insulation between the probe and 
the tube wall. More important, since teflon is not wetted 
by water, there is no chance of a short-circuit between the 
exposed end of the wire and the tube wall. The drying time 

of the teflon after passage of a slug of liquid is extremely 

-3 

short, less than 10 seconds as observed with an oscillo- 
scope 

2.2.3 Probe Sensitivity 

A schematic diagram of the probe circuit is shown in 
Fig. 3a. It can be shown that the voltage reading, Vj, is 
related to the resistance of the water-steam mixture in the 
test section, R^, by the equation 

Vi 

Furthermore, It can be shown that the maximum sensitivity, 
defined by a maximum in Is obtained when 

J- = J^+ J- 

R,. R„ Rt 



( 2 ) 
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At this maximutn sensitivity, the instrument voltage, is 

related to the battery voltage, V^, by the equation 



V = — 

''12 



(3) 



The dependence of V^- on R , for V-o = 6 volts and = 

-L VJ JD 1 

10 ohms, is indicated in Fig. 3b. The lowest value of R^ 

occurs when saturated liquid is flowing in the test section. 

During the course of the experim.ents „ R for saturated liq- 

w 

c; 5 

uid was found to vary betvjeen 10-^ and 2 x 10 ohms. For the 
lower limit of R^ = 10^ ohms, it can be seen that for R^ = 
lO'^ ohms, the maximum sensitivity will occur in the liquid 
region. For R^ = 00 the maximum sensitivity will occur at a 
much higher R^ corresponding to a flow pattern with a higher 
void fraction. Consequently, by adjusting R^, the maximum 
sensitivity can be chosen to correspond to a flow regime 
boundary. However, since these boundaries are not known in 
advance, it is desirable to maintain R^ at a fixed value for 
each run. It was found most convenient to adjust R^ so that 
maximum sensitivity was obtained in the bubbly region. This 
could not always be achieved due to the restriction imposed 
by the internal resistance of the oscilloscope. This, however 
did not affect the results to any noticeable degree since the 
variable of Interest was the shape rather than the value of Vj 

2 .3 Probe Interpretation 

The output of the probe for typical runs is described 
below. Throughout the interpretation phase of the experi- 
ment, Polaroid photos were simultaneously taken of the exit 
flow and of the oscilloscope trace. 
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2.3.1 Velocity Effect on Probe Signal Amplitude 

With stationary water In the test section, a large 

signal was obtained. When the pump was started and water 
began to flow, the signal decreased by one to two volts. 

As the water velocity Increased, the voltage decreased, 
although the decrease was much less than that observed at 
startup. Other experimenters^^^^ also observed this effect 
and speculated that this is an "electrokinetic” effect, 
associated with the charge-carrying capacity of the water. 

2.3.2 Temperature Effect on Probe Signal Amplitude 

When the power was applied to the test section, the 

bulk fluid temperature Increased and the voltage trace in- 
creased. This voltage increased to a maximum as the satu- 
ration temperature was reached. This was expected since 
the conductivity of the water increases with temperature, 
causing a greater voltage drop across the variable resis- 
tor, . To Insure that this effect was only due to temp- 
erature, power was shut off and the water temperature was 
increased by means of the preheaters. The same effect was 
observed . 

2.3.3 Flow Regimes Related to Probe Signal 

The number of distinct flow regimes of two phase flow 
was shown to be somewhat arbitrary. Four main flow 
patterns were recognized in the present experiments: 
bubbly, slug, annular, and wispy annular flow. The physi- 
cal or geometric description of the observed flow patterns 
is very similar to that given in Section 1.2. 
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Bubbly flow is characterized by a continuous liquid 
phase with small vapor bubbles dispersed somewhat uniformly 
throughout the tube . The characteristic bullet-shaped 
bubbles of slug flow could not be identified under the 
present high heat flux conditions* Slug flow was recog- 
nized whenever large ^ irregular bubbles, whose size was on 
the order of the diameter of the tube had somewhat continu- 
ous slugs of water between them. Annular flow consisted of 
a thin film of liquid, usually quite turbulent, flowing along 
the wall, while vapor and fine droplets flowed in the center 
of the tube. This has often been called spray-annular flow, 
but for the purposes of this report, annular will be used. 

The character of the core flow was inferred from the charge 
carrying spray effect which is discussed in Sec, 2.3.4. 

Wispy annular flow was characterized by a continuously in- 
creasing void fraction without the distinct characteristics 
of the above three flows. 

Since the classification of these flow regimes is 
rather qualitative, different individuals may interpret 
the flows differently. The data in Figs. 4, 5, and 6 were 
chosen so as to Indicate the semi-quantitative definitions 
of the flows which were used throughout this study. 

Figure 7 presents a summary of the probe traces for the 
various flow patterns. 

Bubbly Flow In Photo I in Fig. 4, a 2.8-volt signal with 
very small fluctuations is observed » These small fluctua- 
tions are due to bubbles which can be seen in the sight 
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section. Figure 7a shows the schematic of the oscillo- 
scope trace for bubbly flow. Incipient boiling could not 
be detected by the probe due to the small amount of vapor 
generation. It could not be detected with the sight sec- 
tion either, because these bubbles recondensed in the 
subcooled bulk flow before visual observation could be 
made . 

Bubbly-to-Slug Transition The transition from bubbly to 
slug (BTS) flow was first defined when the small pips 
previously seen in bubbly flow became large enough to pro- 
duce a 1 .0 - 1.5-volt deflection. Photo II graphically 
shows a large bubble just leaving the heated section. 

This bubble caused a 1.3-volt deflection on the scope 
trace. This transition, schematically shown in Fig. 7b, 
was well defined and reproducible. 

Slug Flow Slug flow was considered to be present when 
the probe signal alternated between all-liquid and all- 
vapor levels. In Photo III it is noted that clean 
bubbles and slugs, of the variety found in adiabatic flow, 
do not exist . The bubble is like a vapor void without a 
definite head , or tail, and the slugs have small bubbles 
dispersed throughout. At high flow rates, the slug trace 
was somewhat different, as seen in Photo IV. The frequency 
of the full deflection is increased dramatically and the 
flow appears to be very agitated. Figures 7c and 7ci show 
the first full slug and slug flow, respectively. 

Slug-to -Annular Transition The transition criterion for 



the slug to annular (STA) was more subjective than the 
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other criteria o From Fig. 7d it is seen that for slug flow, 
the voltage deflections average out between and 0. 

In Fig. it appears that the slugs are thinner and a 
larger percentage of vapor is passing the probe . The STA 
transition was defined when the voltage deflections appeared 
to rise from 0 volts up, rather than from down. An 

example of this flow is seen in Photo V. It is noted that 
it would be very difficult to distinguish this flow from 
annular flow, particularly at high G, using the sight sec- 
tion alone. As test-section power increases further 
(Photo VI, Fig. 7 ^)) the voltage trace is similar to that of 
the initial description of the STA transition except that 
the peak voltages are not full deflections. This flow is 
still the transitional flow from slug to annular. Physi- 
cally, it appears as if most of the tube consists of annular 
flow, but much of the liquid transport is still in the very 
rapidly moving slugs. 

Annular Flow With greater power, the voltage deflections 
become smaller, thus indicating that the annular flow region 
is being approached. Annular flow was defined to be that 
flow which would give an average deflection no greater than 
0.5 volts (Fig. 7g)o This definition did allow larger de- 
flections to occur, but the frequency of these large deflec- 
tions was greatly reduced. At low velocity annular flow 
was well defined since the voltage was practically zero 
(Photo VII). Photo VIII is another example of annular flow. 
The Increased velocity of the vapor core caused large wave- 
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the initial description of the STA transition except that 
the peak voltages are not full deflections. This flow is 
still the transitional flow from slug to annular. Physi- 
cally, it appears as if most of the tube consists of annular 
flow, but much of the liquid transport is still in the very 
rapidly moving slugs. 

Annular Flow With greater power, the voltage deflections 
become smaller, thus Indicating that the annular flow region 
is being approached. Annular flow was defined to be that 
flow which would give an average deflection no greater than 
0.5 volts (Fig. 7g ) 0 This definition did allow larger de- 
flections to occur, but the frequency of these large deflec- 
tions was greatly reduced. At low velocity annular flow 
was well defined since the voltage was practically zero 
(Photo VII). Photo VIII is another example of annular flow. 
The increased velocity of the vapor core caused large wave- 
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like motions on the surface of the annular film. These 
waves appeared as large deflections on the oscilloscope . 

As G was increased, these waves increased in frequency and 
the film had a milky appearance, probably due to the in- 
creased entrainment and deposition of droplets. Very 
little could have been said about the interior flow, l.e. 
waves, if only photography had been used in this study. 

Wispy Annular Under certain conditions, i.e. = 200°F, 

P = 20-30 psla, G =1.6 - 2.0 X 10^ Ibm/hr-ft^, a varla- 
tion of the above flow regimes was noted. As power was 
Increased beyond bubbly flow, a strong bubbly trace 
travelled downward from to 0 volts. Slug flow did 

not appear and in its place was a flow which can be called 
wispy annular^^^ or froth^^^. The behavior of the trace 
is schematically shown in Figs. 7b, 71^ and 7j and in the 
corresponding Photos IX, X, and XI. As the void fraction 
is Increased, the resistance between probe and wall 
becomes greater until annular flow is reached. It would 
be possible to describe many flow patterns in this high- 
velocity transition between bubbly and annular flow. 

Froth flow would probably be a more appropriate descrip- 
tion of the beginning of this transition, while wispy annu- 
lar could describe the latter part of the transition. 

Since the probe indicated that the transition was quite 
homogeneous throughout, it appears that a single classi- 
fication is adequate . 
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2,3.4 Charge Carrying Spray Effect 

Photo XII again shows annular flow, but the oscillo- 
scope trace is radically different in that negative traces 
up to one volt are seen. This type of trace was not ob- 
served in high pressure studies^^^^^^ and initially caused 
quite a bit of trouble in interpretation. With the probe 
circuit battery removed and the oscilloscope connected 
directly across the probe, a fluctuating negative voltage 
was still observed. The magnitude and frequency of this 
fluctuation increased as the flow rate was increased, even 
though no voltage source was present. A possible pickup 
effect on the probe from the voltage across the test tube 
was discounted because the generator output voltage did 
not have such large fluctuations. 

The explanation of this negative voltage effect is 
believed to be as follows. As annular flow is approached, 
small droplets of water are torn away from the liquid 
film and are transported as a spray in the core. By some 
physical effect these droplets become negatively charged 
upon detaching from the bulk flow, thus leaving the film 
positively charged. The droplets travel down the tube and 
rejoin the film and lose their charge, escape through the 
exit and lose their charge in the loop, or deposit their 
charge on the probe. This last possibility appears to be 
the source of the negative voltage observed on the scope. 
Professor D. P. Keily^^^^ has observed such effects with 
aerosol-type flows and indicated that voltages up to 
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several volts could be observed if there were enough drop- 
lets. The mechanism that causes the detached particle to 
have a negative charge has not been identified. It appears 
that this effect is quite important and should receive more 

( IQ ) 

attention. Vonnegut, et al,' indicate that such charged 
water droplets are produced naturally in clouds and thunder- 
storms and as spray from lakes and oceans. 

This negatively-charged droplet theory is substanti- 
ated by several observations. When the water was particularly 
pure, i .e , 3 megohm-cm and above, this effect was quite 
strong. However, after adding salt to the water (30-40 
ppm), the negative voltage was not observed. The sketches 
k and 1 in Pig. 7 show this purity effect. 

The more conductive water allows more of the negatively 
charged particles to transfer their charge to the liquid 
film rather than depositing it on the probe . This also 
explains why the Dynatech workers, who used a similar 
probe, did not observe this effect. They were operating 
at 500 and 1000 psia and had a correspondingly higher temp- 
erature level. The conductivity of water Increases greatly 
as the temperature increases; for example, the conductivity 
of pure water due to ionization is about a 1000 times 
greater at 546*^? than at Consequently, this con- 

ductive water allowed the negative charge to leak into the 
liquid film before hitting the probe <, With this theory, it 
would be expected that the negative voltage would increase 
with Increased mass flow rate , This effect was also obser- 
ved, as stated earlier. 
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Several methods were attempted to eliminate this effect 
electronically by using a higher battery voltage, and by 
using an ac high-frequency voltage source . The first 
method reduced the relative voltage of the effect but did 
not eliminate it. By means of filter circuits and high 
frequency ac, it was hoped to eliminate the spurious nega- 
tive charge and retain only the desired signals. However, 
the filters and bias circuits also removed the desired 
fluctuations. With experience it was found that proper 
interpretations of the flows could be made even though the 
negative voltage was present. 

2.3.5 Reproducibility and Reliability of Probe Output 

About a dozen runs under various conditions were 
repeated for a second time. When these runs were compared 
to the original ones, it was seen that they were almost 
identical. It was concluded that any run could be repeated 
with the expectancy of getting very similar results. 

The sight section observations were quite helpful in 
Interpreting the probe traces. However, when the bulk 
temperature was moderately subcooled, the vapor condensed 
rapidly in the sight section and it was necessary to look 
into the exit of the heated section to get an accurate 
visual indication of the flow pattern. When the void frac- 
tion was large, the flow was usually so disturbed that it 
was difficult to visualize the flow pattern. In both of 
these situations the probe was able to distinguish the 
important features of the flow. 
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2 .4 Test Program 

After the probe calibration and interpretation runs 
were completed, the actual flow-regime investigations were 
started. The following table lists all the sets of runs 
completed in this experimental program. 

All the data are presented on flow regime maps. The 
coordinates of mass velocity (g) and mass quality (x) were 
chosen because it appeared that they were best suited in 
presenting the various effects studied. In addition, they 
are usually the most convenient for engineering applica- 
tions. Researchers in adiabatic two-phase flow frequently 
make use of volume flow rate and volume quality coordinates 
They could not be used for this study since nonequilibrium 
effects dominated and these quantities cannot be calculated 
Several runs are also presented in terms of heat flux and 
quality to demonstrate a particular effect. 

IXie to the manner in which the loop operates, it was 
difficult to control the exit pressure very accurately. 
Consequently, pressure comparisons were made at low 
pressures where the exit pressure varied from 20-35 psla 
and at high pressures between 90-100 psia , These small 
pressure differences did not affect the results to any 
appreciable degree, in comparison to the effects due to the 
large difference between high and low pressures. 

2 .5 Data Reduction 

All data were reduced by machine computation. The IBM 
709 ^ Computer at the MIT Computation Center was used. The 
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program is described in Appendix C. 
determined by the formula 
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The heat flux was 



( 4 ) 



The equilibrium exit quality is given by 
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q/w + h^ - h^ 
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Chapter III 

PRESENTATION AND DISCUSSION OF RESULTS 
3 .1 Effects of Variables on F1o\m Regimes 

Experiments were conducted to investigate the effects of 
velocity., inlet temperature, exit pressure, quality, purity, 
length, diameter and orientation on the flow regime boundar- 
ies, Figures 8 through 20 present all of the flow regime 
data in terms of mass velocity versus quality. In accordance 
with the discussion in Chapter II, the following flow patterns 
were Identified: bubbly, bubbly-to-slug transition, slug, 

slug-to-annular transition, annular, and wispy annular. 

It is seen in all the maps that the bubbly and slug 
flows occurred primarily in the subcooled region, while full 
annular flow was developed at 2 to 4 percent quality. The 
two-phase flow is therefore in a highly nonequilibrium state 
and a considerable volume of vapor coexists vjith the sub- 
cooled liquid. 

Prom all the maps, it is observed that the transitions 
occur at greater subcooling or lower quality as the velocity 
of the flow is increased. 

The influence of the remaining variables on the flow 
regime boundaries can best be illustrated by composite plots. 
At this point it is desirable to consider only the major 
flow regimes. It is noted that the transition regions 
occasionally cover a broad enough range of quality so that 
they could possibly be classified as flow regimes. However, 
practical considerations dictate that the number of regimes 
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be kept to a minimum. Boundary lines were drawn through 
the center of the transition regions and these lines were 
then transferred to the appropriate composite plots. These 
plots are presented in Figs. 21 through 2?. 

Figures 21 and 22 Indicate that the transitions occur 
at higher quality as the inlet temperature is increased 
from 100 to 200°F. The transition bands for T^ = 200°F 
are much narrower than those for T^ = lOO^F (Figs. 10 and 
12) . In addition, it is noted that slug flow was not ob- 

P 

served above G = 1.2 x 10 Ibm/hr-ft with the higher inlet 
temperature. These latter effects were also observed for 
the smaller diameter test section (Fig, 19), but were not 
observed at higher pressure (Pig. l4). 

At T^ = lOO^P, at 30 and 100 psia, particularly at 
high G, the transition from slug to annular flow was notice- 
ably more violent causing the test section, the connecting 
ih:ibber hoses, and piping to vibrate considerably. While 
testing the L/D = 30 tube, violent vibrations were also 
observed at this transition. This tube burned out in this 
region in spite of the fact that the heat flux was not 
near the critical flux reported by Lopina,^^^^ 

Pressure had a considerable effect on the location of 
the flow regime boundaries. Increasing the pressure from 
30 psia to 100 psia shifted the flow regimes to substan- 
tially lower quality. Figure 23 shows this pressure effect 
at = 100 and 200°P. 

In order to eliminate the negative voltage problem 
previously described, salt was added to increase the conduc- 
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tivity of the water. The effect of this dissolved impurity 
on the flo^^^ regimes is illustrated in Pig, 24, The BTS 
transition for very pure water (2 ,5 megohm-cm) occurs at 
slightly higher quality than the transition for the 1 and 
9 PPM water. There was no discernible effect in the STA 
transition. Conductivity variations between 1 and 3 
megohm«cm produced no noticeable shift in the flow regimes, 
Length-to-diameter ratios of 30, 60 and 90 were 
examined. From Fig, 25, it is seen that the flow regimes 
are quite similar, although for an L/D = 30 the transitions 
are at slightly lower quality. 

For the same T^ there was no discernable difference in 
the flow regime maps for the two diameters (Fig, 26), 

In Pig, 27 it is seen that the transition lines have 
the same general shape for the horizontal and vertical test 

sections. There is a slight shift to lower quality for the 

6 2 

vertical section. At velocities below 0,4 x 10 Ibro/hr-ft , 
slight stratification was visually observed for the bubbly 
and BTS transitions. Since the vertical test section 
showed only a minor flow regime variation from the horizon- 
tal section, it is evident that this stratification did not 
affect the probe output at the particular position, i,e, 
center, that the probe was in, A modification of the 
probe design, whereby the probe could be moved radially 
across the tube during operation, would allow examination 
of the flow very close to the walls and thus would detect 
any stratification. 
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3 .2 Explanation of Results 

The various effects noted above appear to be very 
closely related to Incipient boiling, heat flux, and tur- 
bulence of the flow. To explain this relationship, the 
requirements for the transition to slug flow are first 
examined « 

C 21 ^ 

Radovich and Moissis' ' have shown that the frequency 
of bubble collisions is very high in two-phase bubbly flow. 

A relatively small number of these collisions results in 
coalescence. As more bubbles coalesce, the bubbles 
increase in size to such an extent that slug flow is obser- 
ved. The time that it takes to cause the first identifiable 
existence of slug flow has been called the agglomeration 
time. The agglomeration time is thus determined by the 
collision frequency and the probability of coalescence per 
collision. The latter quantity was found to be strongly 
dependent on the purity of the liquid. Since the liquid 
used in this study was maintained at relatively constant 
purity for most of the experiments, the agglomeration time 
depended primarily on the collision frequency. A decrease 
in the agglomeration time would mean a shifting of the BTS 
transition to lower quality. This could be accomplished 
by increasing the number of bubbles and by causing greater 
turbulence in the flow. 

The requirements for incipient nucleation in a forced- 

convection surface-boiling system have been considered in 

( 22 ) 

detail by Bergles and Rohsenow,' When their correlation 
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was applied to the present conditions (Appendix D), it was 
found that incipient boiling occurred at greater subcooling 
as G increased, decreased, and increased. Further- 
more, it can be seen from Figs. 28 to 31 that the heat flux 
level increases as G increases, decreases,, and 
increases. As the heat flux increases there is an increase 
in the bubble population. Bowring^^^^ has concluded that 
as heat flux and velocity increase, there is an increase 
in the subcooling at which detatchment of the bubbles from 
the heated surface occurs. This is due primarily to a 
general increase in turbulent mixing under these conditions. 
This enhanced mixing should promote the agglomeration of 
bubbles . 

Thus, on the basis of incipient boiling, heat flux, 
and mixing considerations, one would expect a shift in the 
BTS transition to greater subcooling as G increases, T^ 
decreases, and P increases. These trends are observed in 
the present data, as was discussed in the preceding section. 
Since the heat flux level Increases as the tube length 
decreases, the shifting of the BTS transition to greater sub- 
cooling as length is decreased is also accounted for. 

Regarding the other transitions, it is reasonable to 
expect that if the BTS transition occurs at lower quality, 
the STA transition should also occur at lower quality. In 
particular, the increased bubble agitation at the higher 
heat fluxes would be expected to break up the slugs, thereby 
promoting annular flow. 
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The limited occurrance of wispy annular flow can be 
related to the above arguments. It is recalled that this 
regime occurred only for high inlet temperature and low 
pressure. Under these conditions the heat flux level is 
low and there are fewer bubbles and also less tendency for 
the bubbles to agglomerate into slugs. The core flow 
remains bubbly or forthy in appearance until the void 
fraction is so high that the spray annular configuration 
is reached. In general the wispy annular regime is only 
observed at high velocity where the short transit times are 
less condusive to formation of slug flow. 

At this point it is worth mentioning that Radovcich 
and Moissls predicted that the BTS transition would be much 
more rapid for pure water. As discussed in the preceding 
section, the present tests indicated that this transition 
for the demineralized water occurred at about the same 
quality as that for the water with dissolved salt. 

Their observation was for adiabatic flow; thus, with heat 
addition the other mechanisms appear to overshadow the 
impurity effect . 

3.3 Interpretation of Critical Heat Flux Data by Means of 

Flow Regimes 

The present flow regime observations should be useful 
in interpreting heat transfer and pressure drop data for 
boiling water at low pressures. As an example of the 
application of the present data, one of the important 
effects in bulk boiling burnout will be considered. 
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Recent studies of bulk boiling at low pressure have 
confirmed that the critical condition coincides with the 
disappearance of the liquid film from the heated surface 
The film can be depleted by evaporation and entrainment., or 
it can be disrupted by bubble nucleatlono In addition to 
these local effects It appears that unstable upstream flow 
patterns can produce a disruption of the annular film 
further downstream* 

Lopina^^^^ found that the burnout heat flux increased 
as the inlet temperature increased* Typical data illustra- 
ting this effect are plotted in Fig* 33* The flow regimes 
for these conditions are superimposed on the operating lines* 

First consider operation at high inlet tem.perature as 
shown in Fig* 33* The majority of the tube is in annular 
flow with only moderate increases in heat flux* The slug 
flow region of the tube becomes increasingly violent, but 
also moves toward the inlet, as the heat flux is Increased. 
This unstable slug flow region should have little effect on 
the critical heat flux, which for these flow conditions is 
probably caused by nuc lea tion<- induced breakup of the film* 
Now, as the inlet temperature is decreased, the slug flow 
region remains in closer proximity to the exit of the tube 
and the disturbance created within this region would be 
expectedi to propagate downstream and affect the film at its 
thinnest point, which is at the exit of the tube. A 
significant reduction in the critical heat flux is thus 
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Observed as the inlet temperature is increased. 

With short test sections the operating region will lie 
in the vicinity of zero quality. In this case the unstable 
slug flow regime is at the exit to the test section and 
relatively low burnout fluxes, noted in the present L/D = 
30 tests, for example, would be expected. 
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Chapter IV 
CONCLUSIONS 

4 .1 Conclusions 

The conclusions of this investigation can be summarized 
as follows: 

1) The electrical resistance probe is exceptionally 
well suited for flow regime studies » It is simple, inexpen- 
sive, and provides reliable and consistent results o For 
low pressure studies, where nonequilibrium effects dominate, 
the probe is more reliable than an exit sight section » 

2) The Interesting phenomenon of the spray particles 
carrying a negative charge has been observed. The cause of 
this effect is not yet known, but various researchers are 
investigating the problem. Although it influenced the 
probe output under certain conditions, it did not Inhibit the 
flow regime analysis, 

3) Plow regime observations were made for low pressure 
water over a wide range of flow conditions. The major flow 
regimes observed were bubbly, slug, annular, and wispy annu- 
lar flow. The flow regime maps show that bubbly, transition 
to slug, slug, and transition to annular flow occur in the 
subcooled region under nonequilibrium conditions while a 
fully developed annular flow exists by the time 2-4 percent 
quality is reached. 

4) The observed effects of the parameters on the 
flow regimes are summarized in the table below. 
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Table II 



Parameter 


Direction of 
Change in 
Para-meter 


Direction of Shift of 
Flow Regime Boundaries 


Mass Velocity 


Increase 


towards greater subcooling 


Inlet Temper- 
ature 


Increase 


towards less subcooling 


Pressure 


Increase 


towards greater subcooling 


Purity 


Decreased 


negligible change 


Length 


Increased 


slight change towards less sub 
cooling 


Diameter 


Decreased 


negligible change 


Orientation 


Horizontal to 
vertical 


negligible change 



5) The influence of velocity, temperature, pressure and 
length on the flow regime boundaries can be explained in terms 
of incipient boiling, heat flux, and turbulent mixing consid- 
erations . 

6) The present observations were applied to bulk boil- 
ing burnout data . It was found that the effect of inlet temp- 
erature on burnout can be related to unstable upstream flow 
patterns . 
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appendix A 

Burnout Data 



Date ; 19 and 21 July I965 
= lOO^F 

i.d. = 0.242 in. 

L/D = 60 



Run 

1 

2 

3 

4 

5 

6 

7 

8 





G X 10 "^ 


X 


(q/A)cr X lO” 


(psla) 


(Ibm/hr-ft^) 


percent 


(Btu/hr-ft^) 


27.0 


0.5 


33.0 


0.962 


28.0 


0.4 


37.0 


0.843 


43.5 


0.8 


30.6 


1.51 


34.0 


0.8 


19.0 


1.12 


37.5 


1.2 


11.6 


1.36 


36.0 


1.2 


10.9 


1.32 


51.0 


1.6 


10.76 


1.81 


45.0 


1.6 


8.35 


1.64 


48.0 


2.0 


4.15 


1.71 



9 
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APPENDIX B 

Photographic Procedures 



For the photographic study, simultaneous pictures were 
taken of the oscilloscope trace and of the exit flow in the 
glass sight section. This was accomplished by shooting both 
pictures manually at the same time at the desired flow 
conditions . 

The Polaroid Oscilloscope Camera for the Hewlett-Packard 
Company Model 130B/BR Oscilloscope was used. The film was 
Polaroid ASA 3000 film and the camera settings were f 2.8 
and speed of 1/30 sec . 

The sight section was normally illuminated by a 
General Radio Co, Strobotac . The Strobotac was always off 
when actually taking photos . 

For the sight section photos, Polaroid ASA 200 Film was 
used in order to obtain better quality pictures* A General 
Radio Co. Type 1530-A Microflash was used to illuminate the 
exit test section. The 2-microsecond light pulse effectively 
stopped the motion of the flow in the sight section. Various 
backgrounds for the sight section were tried with black 
velvet giving the best results. The Polaroid camera was 
mounted on a tripod and the +2 and +4 Polaroid close-up 
lenses were used. The camera settings were f.8 and 1/30 sec. 

Below is a sketch of the lighting setup viewed from 
above. The microflash was approximately 12 inches below the 
test section and was inclined so that the light was pointed 
directly at the sight section. 
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APPENDIX C 



Computer Program for Data Reduction 
The experimental data obtained was reduced In accor- 
dance with the equations in Section 2 <.5 by a program written 
in FORTRAN and run on the MIT Computation Center IBM 709^ 
Computer . 

The input data consisted of the run identification 
number, the exit pressure, voltage across the test section, 
shunt voltage (current through test section), enthalpy of 
the incoming fluid, and the mass flow rate. 

The program calculated the rate of heat transfer, heat 
flux, rate of heat transfer per pound mass of water, enthal- 
py of saturated liquid, latent heat, and quality. 

Below is a list of the variables as used in the program 
and their physical counterpart . 



PROGRAM VARIABLE 
I, J. K 
P EXIT 

E 

CUR 

HIN 

W 

A 

HADDED 

HTPLUX 



PHYSICAL COUNTERPART 

Numeral run identification numbers 

P - the pressure at the test section 
^exit (psia) 

V - voltage across test section (volts) 

I - shunt voltage (mv) 

h. - enthalpy of incoming subcooled 
water 

w - mass flow rate (Ib/hr) 

A - heat transfer area (Changed this 
value for each L/D and d) (ft^) 

q - rate of heat flow (Btu/hr) 

q/A - heat flux (Btu/hr-ft^) 
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BTULBM 


q/w - rate of heat flow per Ibm of 
water (Btu/hr-lbm) 


ENFPX 


hfe - enthalpy of saturated liquid at 
exit (Btu/lbm) 


ENFGPX 


h^g - latent heat at exit (Btu/lbm) 


X 


X - quality 



The enthalpy of saturated liquid (h^) as a function 

f 2*5 ) 

of pressure was obtained by Todreas^ by a power series 
fit of the Keenan and Keyes* values. For P < 200 psia 



h.(P) = 2 a.P^ 

^ i=o ^ 

where 

= 1.1322734(10^) 
a^ = 6.3204790(10°) 
ag = -1.4742752(10"^) 
ag = 2.5403593(10"^) 
aj^ = -2.8788220(10"5) 

The latent heat (h^^) wa 
Keenan and Keyes* values on a 



a^ = 2.1239014(10"’^) 
ag = -1.0067598(10"^) 

= 2.9480958(10"^^) 

3g = -4.8430355(10"^^) 

= 3.4075972(10"^®) 

obtained by fitting the 
power series. For P <400 psia, 



h_(P) = 2 a pi 

i=o ^ 

= 9.9704457(10^) 
a^ = -2.1762821(10°) 
ag = 1.9558791(10"^) 

83 = -1.2610018(10"'^) 

= 5.0308804(10""^) 

The values of the heat 



= -1.2499843(10"^) 
ag = 1.9322768(10"^^) 

a^ = -1.8023251(10"i®) 
ag = 9.2680092(10"i^) 
a^ = -2.0148643(10"^^) 
transfer area for the differ- 



ent tubes are: 
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TUBE 


A 


0.242 in. 


i.d. L/D = 30 


0.038274 


0.242 


60 


0.076548 


0.242 


90 


0.114822 


0.094 


60 


0.0154 



The output data consisted of all the input variables 
and the computed quantities. 

The program used for the L/D = 60, D = 0 .094-ln. test 
tube computation follows. 

COMPUTATION PROGRAM 



1 READ 2, I, J, K, PEXIT, E, CUR, HIN, W 

A = .0154 

HADDED = 205.36 * E * CUR 
HTFLUX = HADDED / A 
BTUEBM = HADDED / W 
P = PEXIT 

5 ENPPX = 1.1222734E2 + 6.3204790*P - 1 .4742752E-1*(P»»2) 

1 +2.5403593E-3 *(P**3) -2 .8788220E-5*(P*»4) +2.1239014E-7*(P**5) 

2 -1.0067598e-9*(p**6) +2.9480958e-12*(p**7) 

3 -4.8430355E-15*(P**8) +3 .4075972E-l8*(P**9) 

6 ENFGPX=9 .9704457E2 -2.1762821*P +1 .9558791E-2*(P**2) 

1 -1.2610018e-4*(p**3) +5.0308804e-7*(P**4) -1.2499843E-9*(P**5) 

2 +1. 93227 68E-12*(P**6) -1.8023251E-15*(P**7) 

3 +9.2680092E-19*(P**8) -2.0148643E-22*(P**9) 

X = (HIN + BTULBM - ENFPX)/ENFGPX 

PRINT 3, I, J, K, PEXIT, E, 

I CUR, HIN, W, HADDED, HTFLUX, BTULBM, X, ENFPX, ENFGPX 
GO TO 1 

2 FORMAT (315,5E10.2) 

3 FORMAT (7HORUN NOI3, 1H-, 13, 1H-,13,5X, 

1 6HPEXIT=E15.6,5X,2HE=E15.6,5X,4HCUR=E15.6/5X, 

2 4hhin=E15.6,5X2HW=E15.6,5X, 7 

3 HHADDED=E15 • 6/5X, 71flITFLUX=EI5 .6, 5X, 

4 7HBTULBM=E15 .6,5X,2HX=E15 .6/5X 

5 6HENPPX=E1 5 . 6 , 5X , 7HENFG PX=E15 . 6 ) 

END 
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appendix D 

Incipient Bolling Analysis 

The effects of velocity, temperature, and pressure on 
the Incipient boiling heat flux, (q/A)^, can be examined by 
means of the correlations given in Ref, 22, This correla- 
tion for water is: 

2,30 

T ^,0234 

(q/A)^ = 15.60 (Dl) 

For a uniformly heated tube the heat flux can be calculated 
from: 



(q/A) =h(T^-Tjj) =h (VTg) + (Tg-T^) (D2) 

The saturation temperature is known from the exit pressure 
The bulk temperature is found by the First Laws 

which leads to: 

'^out = . (d4) 

Since the probe is at the exit, only the conditions at the 

exit are of interest and therefore T, = T . . The heat 

b out 

transfer coefficient, h, is found using McAdams* correla- 




where the properties are evaluated at the bulk fluid tempera- 
ture. For the actual numerical analysis, the values of these 
properties came from ElWakil By assuming a value of 








\ 
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(q/A), Tj^ and hence h can be calculated. Then using (D2), 
(T^-T^) is found. The values of (q/A) and (T^-T^), for 
various, G, T^, and can then be found and plotted. 
Equation (Dl) for the various pressures Is then superimposed 
on these plots. The intersections of the (q/A) vs. 
curve and the (q/A)^ curve Is the point at which incipient 
boiling is predicted. 

This analysis was carried out for the conditions that 
were presented in Figs. 28 through 31, i.e. 





G X 10"^ 
Ibm/hr-ft^ 




°F 


^e 

(psla) 


a) 


0.4 


100, 


150 , 200 


30 


b) 


0.4 


100, 


200 


100 


c) 


2.0 


100, 


150, 200 


30 


d) 


2.0 


100, 


200 


100 



In Fig. 32a is a representative plot for the conditions in 
(c). To visually relate the predicted (q/A)^ to (x), opera- 
ting lines obtained from the formula: 



•*- hjn - 



(d6) 



are plotted on coordinates of (q/A) vs x. The (q/A)^ found 
from the analysis are plotted on these operating lines and are 
connected by a smooth line. This line then is the locus of 
predicted (q/A)^ for various inlet temperatures. This plot 
for conditions (c) above is presented in Fig. 32b. Examin- 
ing Fig. 32b it is seen that incipient boiling occurs at 
greater subcoolings for lower T^ . From the analysis of 
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conditions a, b, d, it was also found that incipient boiling 
occurs at greater subcooling for larger pressures and, keeping 
P and constant, at higher G. 



MOTOR - GENERATOR 
POWER SUPPLY 
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FIG. 1 SCHEMATIC LAYOUT OP TEST FACILITY 






BURNOUT TEST SECTION 



TEST SECTION EXIT 
WITH SIGHT SECTION 



FIG. 2 DETAILED VIEV/ OF TEST SECTIONS 
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PIG. 3b CHARACTERISTICS OP PROBE CIRCUIT 



VoltL 










PHOTO I 

P-*'30, G=0.^'Xl0''3 T^=200, x=-0.t'7>-' 
BUBBLY 




PHOTO II 

Pa/30, G=2.0x1u'“, T^=r.OO., x = -2.l1‘. 
BTS TRANSITIOH 








irn^^rruKi 

iiiiiiiii'iiinifiiir/i 

iriumniiiiiifii'i 

IJHUBIIliil 



PHOTO III PHOTO IV 

P'v30, G=0.4xI0’, T^=li;U, X = -. .vi. P'^so, G=2.0xlu' . x = 



JLUG 



oJaUG 



FIG. h COMPARISON OF PROBE AND SIGHT-SECTiu;. ^ . SEKv.'-. ;i^:.S 
(BUBBLY, BUBBLY -TO -SLUG, SLUG) 



HI 
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iHMHiiiHnn 



BSSSiBLl 



0.0''i 



PHOTO V 

P'^30, 0=2.0x10', T^=100, x=0.2::i 
FIRST INDICATION OP STA TRANSITION 



i 




PF.UTO VI 

P'''30, 0=0. ^xlO , l^=10vj, 
STA TRANSITION 



m I T 







HMIHISIMIMIH 


■HMffliUiHiP 








■■■nunni 


mmmisimmmm 
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PHOTO VII 

P'vSO, 0=0.4x10^', T^=200, x=10.99; 
ANNULAR 



PHOTO VIII 
P'-jO, 0=1.2x10', 

ANNULAR 



FIG. 5 COMPARISON OP PROBE AND SIGHT -SECTION OBSERVATIONS 
(SLUG-TO-ANNULAR, ANNULAR) 



Volts 











0 



0 0.04 

^ sec 

PHOTO IX 

P-vSO, G=2.0xl0^, T^=200, x=-0,7o'. 
V/ISPY ANNULAR 




iiiy L 

lllklll'^lliCBHiB 



PHOTO X 

P^30, G=2.0xl0" , T^=200, J. ■■ 
WISPY ANNULAR 




PHOTO XI 

P'^SO, G=2.0xl0^, T^=20U, x=4.6S':. 

WISPY ANNULAR 

FIG. 6 COMPARISON OP PROBE 
(WISPY ANNULAR) 




PHOTO XII 

p^30, G=0.uxl0", V.'.' 

ANNULAR WITH NEGATiVL •/■.lARGE Sr-^HC 

SIGHT -SECTION OBSERVATIONS 
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FIG. 7 SCHEMATIC OF VOLTAGE TRACES FOR VARIOUS FLOV/ REGIMES 
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FIG. d FLOW REGIME MAP - V/ATER PURITY = 9PPM NaCl 
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FIG. 9 FLOW REGIME MAP - WATER PURITY = 0,48 Megohm-cm 
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FIG. 10 FLOW REGIME MAP - V/ATER PURITY = 2.5 Megohm-cir; 
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FIG. 11 FLOW REGIME MAP - T. = 1^0 “F, LOW PRESSURE 
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?IG. 12 PLOW REGIME MAP - T = 200 °F, LOW PRESSURE 
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PIG. 13 PLOW REGIME MAP - HIGH PRESSURE, T. = 100 °P 
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FIG. 14 FLOIV REGIME MAP - HIGH PRESSURE, T. = 200 
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FIG. 15 FLOW REGIME MAP - L/D 
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FIG. 16 FLOW REGIME MAP - L/D = 90 
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FIG. 17 FLOW REGIME MAP - L/D = 90, HIGH PRESSURE 
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PIG. Id FLOV/ REGIME MAP - D = 0.094 in., T. = 100 °F 
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FIG. 19 PLOVJ REGIME MAP - D = 0.09^t in., T. - 200'' 
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FIG. 20 FLOW REGIME MAP - VERTICAL UPFLOW 
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-■IG. .;1 lii-’FECT OF INLET TEMPERATURE ON FLOW REGIME BOUNDARIES - P 50 psia 
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FIG. 22 EFFECT OF INLET TEMPERATURE ON FLOW REGIME BOUNDARIES - 100 psia 
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EC'i’ OF fXIT PRESSURE ON FLOW REGIME oUNDARI 
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FIG. 24 EFFECT OF WATER PURITY ON FLOW REGIME BOUNDARIES 
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FIG. 2^ EFFECT OF LENGTH ON FLOW REGIME BOUNDARIES 
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FIG. 26 EFFECT OF DIAMETER ON FLOW REGIME BOUNDARIES 
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FIG. 27 EFFECT OF ORIENTATION ON FLOV/ REGIME BOUNDARIES 
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FIG. 28 FLOW REGIME MAP ON HEAT FLUX vs. QUALITY CO-ORDINATES 
G = 0.4 X 10“ lbm/hr-ft2, LOW PRESSURE 
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FIG. 29 FLOW REGIME MAP ON HEAT FLUX vs. QUALITY CO-ORDINATES 
G = 0.4 X 106 lbm/hr-ft2, HIGH PRESSURE 
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FIG. 31 FLOW REGIME MAP ON HEAT FLUX vs. QUALITY CO-ORDINATES 
G = 2.0 X 10° lbm/hr-ft2, HIGH PRESSURE 



(q/A) X 10’^ BTU/hr-ft 
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FIG. 32a INCIPIENT BOILING CALCULATIONS 




X - percent 

FIG. 32b SAMPLE OPERATING LINES SHOV/ING LOCUS OF INCIPIENT BOILING POINT 
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